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Abstract: We develop a simple semi-analytic model of the star formation rate (SFR) 
as a function of time. We estimate the SFR for a wide range of values of the cosmo- 
logical constant, spatial curvature, and primordial density contrast. Our model can 
predict such parameters in the multiverse, if the underlying theory landscape and the 
cosmological measure are known. 



Contents 



Introduction 

Model 

U Goals 



2.2 Geometry and initial conditions 



2~3| Linear perturbation growth and halo formation 



2~Q Cooling and galaxy formation 



2.5 Star formation 



Results |T3 

~~J] Our universe |13 

1| Varying single parameters |15 

11 Varying multiple parameters E20 



1. Introduction 

In this paper, we study cosmological models that differ from our own in the value of one 
or more of the following parameters: The cosmological constant, the spatial curvature, 
and the strength of primordial density perturbations. We will ask how these parameters 
affect the rate of star formation (SFR), i.e., the stellar mass produced per unit time, 
as a function of cosmological time. 

There are two reasons why one might model the SFR. One, which is not our purpose 
here, might be to fit cosmological or astrophysical parameters, adjusting them until the 
model matches the observed SFR. Our goal is different: We would like to explain the 
observed values of parameters, at least at an order-of-magnitude level. For this purpose, 
we ask what the SFR would look like if some cosmological parameters took different 
values, by amounts far larger than their observational uncertainty. 

This question is reasonable if the observed universe is only a small part of a multi- 
verse, in which such parameters can vary. This situation arises when there are multiple 
long-lived vacua, such as in the landscape of string theory. It has the potential to ex- 
plain the smallness of the observed cosmological constant [1], as well as other unnatural 
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coincidences in the Standard Model and in Standard Cosmology. Moreover, a model 
of star formation can be used to test the landscape, by allowing us to estimate the 
number of observers that will observe particular values of parameters. If only a very 
small fraction of observers measure the values we have observed, the theory is ruled 
out. 

The SFR can be related to observers in a number of ways. It can serve as a starting 
point for computing the rate of entropy production in the universe, a well-defined 
quantity that appears to correlate well with the formation of complex structures such 
as observers [2]. (In our own universe, most entropy is produced by dust heated by 
starlight [3].) Alternatively, one can posit that a fixed (presumably very small) average 
number of observers arises per star, perhaps after a suitable delay time of order billions 
of years. 

To predict a probability distribution for some parameter, this anthropic weighting 
must be combined with the statistical distribution of parameters in the multiverse. 
Computing this distribution requires at least statistical knowledge of the theory land- 
scape, as well as understanding the measure that regulates the infinities arising in the 
eternally inflating multiverse. Here, we consider none of these additional questions and 
focus only on the star formation rate. 

We develop our model in Sec. |2|. In Sec. |3] we display examples of star formation 
rates we have computed for universes with different values of curvature, cosmologi- 
cal constant, and/or primordial density contrast. We identify important qualitative 
features that emerge in exotic regimes. 



2. Model 
2.1 Goals 

Our goal is not a detailed analytical fit to simulations or observations of the SFR in 
our universe. Superior models are already available for this purpose. For example, 
the analytic model of Hernquist and Springel [4] (HS) contains some free parameters, 
tuned to yield good agreement with simulations and data. However, the HS model only 
allows for moderate variation in cosmological parameters. For example, HS estimated 
the SFR in cosmologies with different primordial density contrast, Q, but only by 10%. 

Here we are interested in much larger variations, by orders of magnitude, and 
not only in Q but also in the cosmological constant, A, and in spatial curvature. A 
straightforward extrapolation of the HS model would lead to unphysical predictions in 



-2- 



some parameter ranges. For example, a larger cosmological constant would actually 
enhance the star formation rate at all times [5]. 1 

In fact, important assumptions and approximations in the HS model become in- 
valid under large variations of parameters. HS assume a fixed ratio between virial and 
background density. But in fact, the overdensity of freshly formed haloes drops from 
about 200 to about 16 after the cosmological constant dominates [6]. Moreover, HS 
apply a "fresh-halo approximation" , ascribing a virial density and star formation rate 
to each halo as if it had just been formed. However, large curvature, or a large cosmo- 
logical constant, will disrupt the merging of haloes. One is left with an old population 
of haloes whose density was set at a much earlier time, and whose star formation may 
have ceased due to feedback or lack of baryons. Finally, HS neglect Compton cool- 
ing of virialized haloes against the cosmic microwave background; this effect becomes 
important if structure forms earlier, for example because of increased Q [7]. 

All these effects are negligible at positive redshift in our universe; they only begin 
to be important near the present era. Hence, there was no compelling reason to include 
them in conventional, phenomenologically motivated star formation models. But for 
different values of cosmological parameters, or when we study the entire evolution of 
the universe, these effects become important and must be modelled. 

On the other hand, some aspects of star formation received careful attention in 
the HS model in order to optimize agreement with simulations. Such refinements 
are unnecessary when we work with huge parameter ranges, where uncertainties are 
large in any case. At best, our goal can be to gain a basic understanding of the 
quantitative behavior of the SFR as parameters are varied. We will mainly be interested 
in the height, time, and width of the peak of the SFR depend on parameters. In 
currently favored models of the multiverse [2,8], more detailed aspects of the curve 
do not play an important role. Thus, our model will be relatively crude, even as we 
include qualitatively new phenomena. 

Many extensions of our model are possible. The most obvious is to allow additional 
parameters to vary. Another is to refine the treatment of star formation itself. Among 
the many approximations and simplifications we make, the most significant may be 
our neglect of feedback effects from newly formed stars and from active galactic nuclei. 
In addition, there are two regimes in our model (high Q and negative A near the Big 
Crunch) where we use physical arguments to cut off star formation by hand. These are 
regimes where the star formation rate is so fast that neglected timescales, such as halo 
merger rates and star lifetimes, become important. An extended model which accounts 
for these timescales may avoid the need for an explicit cutoff. 

1 This effect, while probably unphysical, did not significantly affect the probability distributions 
computed by Cline et al. [5]. 



- 3 - 



2.2 Geometry and initial conditions 

Our star formation model can be applied to open, flat, and closed Friedman- Robertson- 
Walker (FRW) universes. However, we will consider only open geometries here, with 
metric 

ds 2 = -dt 2 + a(tf (dx 2 + sinh 2 X dtt 2 ) . (2.1) 

Pocket universes produced by the decay of a false vacuum are open FRW universes [9]; 
thus, this is the case most relevant to a landscape of metastable vacua. The above 
metric includes the flat FRW case in the limit where the spatial curvature radius a is 
much larger than the horizon: 

a > H- 1 , (2.2) 

where H = a/a is the Hubble parameter. This can be satisfied at arbitrarily late times 
t, given sufficiently many e-foldings of slow roll inflation. 

The scale factor, a(t), can be computed by integrating the Friedmann equations, 

2 1 8ttG N 

H - - 2 - — p, (2.3) 
H + H 2 = -^(p + 3p), (2.4) 

beginning at (matter-radiation) equality, the time t cq at which matter and radiation 
have equal density: 

7T 2 

Pm(tcq) = Pr(t cq ) = ^J^cq > ( 2 -5) 

where 

accounts for three species of massless neutrinos. The temperature at equality is [6] 

T eq = 0.22£ = 0.82 eV , (2.7) 
and £ is the matter mass per photon. The time of equality is 

t eq = . « 0.128 G N 1/2 T- 2 = 4.9 x 10 4 yrs . (2.8) 

The total density and pressure are given by 

p = p r + p m + p A , (2.9) 
P = j-PA- (2.10) 
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Among the set of initial data specified at equality, we will treat three elements as 
variable parameters, which define an ensemble of universes. The first two parameters 
are the densities associated with curvature, G^ 1 a(t eq )~ 2 , and with the vacuum, = 
A/87rGN, at equality. We will only consider values that are small compared to the 
radiation and matter density at equality, Eq. (|2.5| ); otherwise, there would be no matter 
dominated era and no star formation. The third parameter is the primordial density 
contrast, Q, discussed in more detail in the next section. 

We find it intuitive to trade the curvature parameter a(t eq ) for an equivalent pa- 
rameter AN, defined as 

a(t cq ) 



eq ; 



(2.11] 



AN = log 
where 

aminW = ( _ 3 4 ^; 4 „ ) V3 (fir)" 172 = 1.13 x 10 24 m = 37Mpc , (2.12) 

is the experimental lower bound on the curvature radius of our universe at matter- 
radiation equality. (Consistency with our earlier assumption of an open universe re- 
quires that we use the upper bound on negative spatial curvature, i.e., on positive values 
of tt k . We are using the 68% limit from WMAP5+BAO+SN [10], tt k = -0.0050^;^. 
We use the best fit values H = 68.8 km/s/Mpc and Q m = 0.282 for the current Hubble 
parameter and matter density fraction.) 

AN can be interpreted as the number of e-folds of inflation in excess (or short) 
of the number required to make our universe (here assumed open) just flat enough 
(O™ = 0.0011) for its spatial curvature to have escaped detection. Requiring that 
curvature be negligible at equality constrains AN to the range AN > —5. Values near 
this cutoff are of no physical interest in any case: Because too much curvature disrupts 
galaxy formation, we will find that star formation already becomes totally ineffective 
below some larger value of AN. 

We hold fixed all other parameters, including Standard Model parameters, the 
baryon to dark matter ratio, and the temperature at equality. We leave the variation 
of such parameters to future work. 



2.3 Linear perturbation growth and halo formation 

Cosmological perturbations are usually specified by a time-dependent power spectrum, 
V, which is a function of the wavenumber of the perturbation, k. The r.m.s. fluctuation 
amplitude, a, within a sphere of radius R, is defined by smoothing the power spectrum 
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with respect to an appropriate window function: 

g 2 = JLjf° ^"Ij'^a . (2.13) 

The radius i? can be exchanged for the mass M of the perturbation, using M = 
47rp m i? 3 /3. Starting at t cq , once all relevant modes are in the horizon, the time devel- 
opment of a can be computed using linear perturbation theory. In the linear theory, all 
modes grow at the same rate and so the time dependence of er is a simple multiplicative 
factor: 

a(M,t) = Qs(M)G(t). (2.14) 

Here Q is the amplitude of primordial density perturbations, which in our universe is 
observed to be Q ~ 2 x 10~ 5 . G(t) is the linear growth function, and is found by 
solving 

d 2 G _„dG 



dfl +2H— = 4irG N p m G (2.15) 

numerically, with the initial conditions G = 5/2 and G = 3H/2 at t — t cq . This 
normalization is chosen so that G ~ 1 + | a ,^\ near t = t ea , which is the exact solution 

2 a(t eq ) "1' 

for a flat universe consisting only of matter and radiation. For s(M) we use the fitting 
formula provided in [6]: 

s(M) = [(9.1^ 2 / 3 ) _0 - 2 V (50.51og 10 (834 + /i- 1 / 3 ) -92)"°- 27 ] _1/0 ' 27 (2.16) 
with \x = Mj M eq , where 

M eq = G N 3/2 r 2 = 1.18 x 1O 17 M (2.17) 

is roughly the mass contained inside the horizon at equality, and £ = 3.7 eV is the matter 
mass per photon. This fitting formula for s(M) and our normalization convention for 
G(t) both make use of our assumption that curvature and pa are negligible at t eq . 

We use the extended Press-Schechter (EPS) formalism [11] to estimate the halo 
mass function and the distribution of halo ages. In this formalism the fraction, F, of 
matter residing in collapsed regions with mass less than M at time t is given by 

F(< M, t) = Erf ( —J^ ] (2.18) 

with 5 C = 1.68 the critical fluctuation for collapse. 5 C is weakly dependent on cosmolog- 
ical parameters, but the variation is small enough to ignore within our approximations. 

Given a halo of mass M at time t, one can ask for the probability, P, that this 
halo virialized before time t vir . As in Ref. [11], we define the virialization time of the 
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halo as the earliest time when its largest ancestor halo had a mass greater than M/2. 
To find P, we first define the function 

/3(Mi, h, M 2 , t 2 ) = Erfc ( = ( - -^—) ) ■ (2.19) 

Then the desired probability is given by 

P(<t vir ,M,t) = / ^-{M lttyil ,M,t)dM 1 . (2.20) 

iM/2 Ml "Mi 

The virialization time further determines the density, p vir , of the halo. We will 
compute the virial density by considering the spherical top-hat collapse of an overdense 
region of space. Birkhoff 's theorem says that we may consider this region to be part 
of an independent FRW universe with a cosmological constant identical to that of the 
background universe. As such, it evolves according to 

2 _ 87rG N Pm , 8nG N p A k _ 87rG N |p A | ( 1 n k 



„ 2 = ^ + ^_„ = ^ r ^_ ±1 __^ (221) 

/ I I \ 1/3 a 2 

where A = ('-^) and k = 5 „ 3 , — , f— ) k. We neglect the radiation term, which 

\ pm J 8ttG n |pa| V a ) to ' 

will be negligible for virialization well after matter-radiation equality. 

Assuming that k and p A are such that the perturbation will eventually collapse, 
the maximum value of A is obtained by solving the equation H = 0: 

±A 3 max - nA max + 1 = 0. (2.22) 

The time of collapse, which we identify with the halo virialization time t v ; r , is twice the 
time to reach A max and is given by 



*vir = 2 A / _^",_ , J 7 =^ (2.23) 



dA 

87tG n |p a | Jo ^±a 2 -k 



To obtain the virial density, we follow the usual rule that p vir is 8 times the matter 
density at turnaround. In other words, 

= 8-^-. (2.24) 

Note that according to this prescription, p v - u {t) has no explicit dependence on the 
curvature of the background universe. Changing p A results in a simple scaling behavior: 
Pvir (t/t\)/ 1 Pa I is independent of p A , where 

t A = (3/87rG N |pA|) 1/2 • (2.25) 
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Pvir/KI 




Figure 1: The virial density is independent of curvature and has a simple scaling dependence 
on pa- For positive p\, p vlv asymptotes to 16/ja at late times. For negative p\, the virial 
density formally goes to zero at t = nt/^, which the latest possible time the universe can crunch. 
As we discuss in the text, our approximation (and, we shall assume, structure formation) fails 
when the virial density drops below the background matter density, which diverges near the 
crunch. 



In practice this means that one only has to compute p v i r (t) for two values of pa (one 
positive and one negative). In Fig. [I], we show the pA-independent part of the virial 
density. 

The virial density determines the gravitational timescale, t gra v, which we define to 

be 

^grav — r-^f • (2.26) 

V CJNPvir 

This is the timescale for dynamical processes within the halo, including star formation. 
2.4 Cooling and galaxy formation 

Halos typically virialize at a high temperature, T vir . Before stars can form, the baryons 
need to cool to lower temperatures and condense. Cooling proceeds by bremsstrahlung, 
molecular or atomic line cooling, and Compton cooling via scattering against the CMB. 
For bremsstrahlung and line cooling, the timescale for the baryons to cool, t coo \, is 



3T vir mi 

^cool 



(2.27) 



2/ b pX 2 A(T vir )p vir 

with p = 0.6 m p the average molecular weight (assuming full ionization), X = .76 the 
hydrogen mass fraction, and A(T v ; r ) the cooling function, which encodes the efficiency 
of the bremsstrahlung and line cooling mechanisms as a function of temperature. For 
T < 10 4 K, A(T) is effectively zero and cooling does not occur. When T > 10 7 K, 
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bremsstrahlung cooling is the dominant mechanism and A(T) oc T 1 / 2 . For intermediate 
values of T, A(T) has a complicated behavior depending on the metallicity of the 
gas [12]. For our purposes, however, it suffices to approximate the cooling function as 
constant in this range. To summarize: 

for T < 10 4 K, 

A(T) = { A„. for 10 4 K < T < 10 7 Is. ;:2.:Jn ; 

A o(lon<) 1/2 > forl0 7 K<T, 

with A = 10 _23 ergcm 3 s _1 . 

We need to include the effects of Compton cooling separately. In general, Compton 
cooling depends on the temperature of the CMB, T C mb, the temperature of the baryonic 
matter, and the ionization fraction of the matter [13]. But the dependence on the gas 
temperature drops out when T vir ^> Tcmb (which is always the case), and we will 
assume that the gas transitions from neutral to completely ionized at T vir = 10 4 K. 
With these approximations, Compton cooling proceeds along the timescale 

45m e , . 

^comp = A o„ (T VP (2.29) 
47T z (Tt(Jcmb) 



where <jt = ( — ) is the Thompson cross-section. 

The Compton timescale is independent of all of the properties of the gas (provided 
T vir > 10 4 K). So this cooling mechanism can still be efficient in regimes where other 
cooling mechanisms fail. In particular, Compton cooling will be extremely effective in 
the early universe while the CMB is still hot. We define r com p as the time when t C om P is 
equal to the age of the universe. r comp is a reasonable convention for the time beyond 
which Compton cooling is no longer efficient. While we compute r comp numerically in 
practice, it is helpful to have an analytic formula for a special case. In a flat, Pa = 
universe we obtain 

<■-» = (d^r f * 0763 &r g - - °' 9Gyr - (2 - 3o) 

Turning on strong curvature or p\ will tend to lower r comp . 

Before r com p, all structures with T v ; r > 10 4 K can cool efficiently and rapidly, so 
star formation will proceed on the timescale t grav . After r comp , some halos will have 
^cooi < tgrav and hence cool efficiently, and some will have t coo \ > t grav and be "cooling- 
limited" . We take it as a phenomenological fact that the cooling-limited halos do not 
form stars. (One could go further and define a "cooling radius" to take into account 
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the radial dependence of the baryonic density [13], but this effect is negligible within 
our approximations.) 

In the absence of Compton cooling, then, we will consider only halos for which 
^cooi < igrav, which have masses in a finite range: 

M min (t vir ) < M < M max (t vir ) . (2.31) 

The minimum halo mass at time t vir corresponds to the minimum temperature T vir = 
10 4 K: 

%r , , /5xl0 4 K\ 3/ V 3 \ 1/2 , nnnS 

m ^ = {-^]t) ■ (2 - 32) 

where we have used the relation between virial temperature, density, and mass for a 
uniform spherical halo derived from the virial theorem: 

T vir = ^ (yPvi^y 73 . (2.33) 

(A non-uniform density profile will result in a change to the numerical factor on the 
RHS of this equation, but such a change makes little difference to the result.) 

The maximum mass haloes at time t vir satisfy t coo i = t gTay . For T < 10 7 K, this 
mass is given by 

Note that halos with mass M max (t) have virial temperatures less than 10 7 K for t > 
0.24 Gyr, which is well before r comp . This means that in our calculations we can safely 
ignore the decreased cooling efficiency for extremely hot halos. 

2.5 Star formation 

For halos satisfying the criterion t cool < t grav , the hot baryonic material is allowed to 
cool. Gravitational collapse proceeds on the longer timescale t grav . We will assume that 
the subsequent formation of stars is also governed by this timescale. Thus, we write 
the star formation rate for a single halo of mass M as 

dM sinsle Mk fuM 

— % (M, tvir) = A-^— = A-^L- , (2.35) 

0,1 f-grav l^vir ) ^grav l^vir ) 

where Mb is the total mass of baryons contained in the halo and /b is the baryon mas 
fraction. We fix / b to the observed value f h 1/6. 2 



2 It would be interesting to consider variations of /b in future work; see Ref. [14] for an environmental 
explanation of the observed value. 
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The order-one factor A parametrizes the detailed physics of gravitational collapse 
and star formation. The free-fall time from the virial radius to a point mass M is 
(37r/32) 1 / 2 t grav , and so we choose to set A = (32/37r) 1//2 . The final SFR calculated 
for our universe is not very sensitive to this choice. Notice that the single halo star 
formation rate depends only on the mass and virialization time of the halo. 

The next step in computing the complete SFR is to average over the virialization 
times of the halos in existence at time t. Define — j. — (M,t) to be the star formation 
rate of all halos that have mass M at time t, averaged over possible virialization times 
£vi r . Now we need to consider the possible virialization times. 

It may happen that some halos in existence at time t are so old that they have 
already converted all of their baryonic mass into stars. These halos should not be 
included in the calculation. Furthermore, feedback effects reduce the maximum amount 
of baryonic matter that can ever go into stars to some fraction / of the total. (/ is a 
free parameter of the model, and we shall fix its value by normalizing the integrated 
SFR for our universe; see below.) We should then only consider halos with t vir > i min , 
where t min is a function of t only and satisfies 

/ 1 dMf R ^ e \ -1 _ / 

t ^min I n j^j- ^ \M ) *minj I ^^grav(^min) • (2.36) 

Though an explicit expression for t min can be found in some special cases, in practice 
we will always solve for it numerically. 

We still need to account for the fact that some halos existing at t formed at t v ; r with 
masses outside the range M mm (tvir) < M < A / / max (t vir ). Since M min is an increasing 
function of time while M max is a decreasing function of time, this condition on the mass 
range is equivalent to an upper bound on t vir : t vir < t max , where t max satisfies either 

M = M max (t max ) (2.37) 

or 

M = M min (t max ) , (2.38) 

or else t max = t if no solution exists. Recall that the condition M < M max is a condition 
on cooling failure, and so is only applicable for i vir > r comp . Halos with t vir < r comp 
have no upper mass limit. Like t m i n , t max is found numerically. Unlike t min , t max is a 
function of both mass and time. 

Now we can use the distribution on halo formation times deduced from Eq. [2.20| 
to get 

JM av g r 4 f, M f)P 1 

dt vlI . (2.39) 



d * -{M,t) 



dt 



Af h M dP . „„ ' 

t-grav ^vir ) C'^vir 
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The final step is to multiply this average by the number of halos of mass M existing 
at time t and add up the contributions from all values of M. Equivalently, we can 
multiply the normalized average star formation rate, ^ — -£ — , by the total amount of 
mass in halos of mass M at time t and then add up the contributions from the different 
mass scales. The Press-Schechter formalism tells us that the fraction of the total mass 
contained in halos of masses between M and M + dM is j^dM, where F(M, t) is the 
Press-Schechter mass function defined in Eq. [2.181 So restricting ourselves to a unit 
of comoving volume with density po (sometimes called the reference density), we find 
that the SFR is given by 

dM / dt vil —(M,t)- ; b (t vil ,M,t) . (2.40) 

Normalizing the SFR Our model contains the free parameter / representing the 
maximum fraction of the total baryonic mass of a halo which can go into stars. We 
will impose a normalization condition on the SFR for our universe to fix its value. 
Nagamine et al. [15] normalize their "fossil" SFR so that 

t0 p,dt = ^, (2.41) 
o iU 

with t Q = 13.7 Gyr the age of the universe. In other words, the total amount of 
mass going into stars should be one tenth of the total baryonic mass in the universe. 
Due to the phenomenon of recycling, this number is more than the actual amount of 
matter presently found in stars. We find that with / = 0.35, our model satisfies this 
normalization condition. 

Early and late time cutoffs Even after t = t eq , baryonic matter remains coupled to 
photons until t = t rEC ~ 3.7x 10 5 yr. So even though dark matter halos can and do begin 
to collapse before recombination, stars will obviously be unable to form. After t = t Tec , 
the baryons are released and will fall into the gravitational wells of the pre-formed dark 
matter halos. Some time At (of order t Tec ) later, the baryonic matter will have fallen 
into the dark matter halos and can form stars normally. In our model, we account for 
this effect by placing further conditions on £ min . First, if t < 2t Tec then the SFR is set 
to zero. Then, provided t > 2t rec , we compute t m[n according to Eq. |2.36|. If we find 



< 2t rcc < t, then the computation of t\ is flawed and we set t min = t cq to 
reflect the fact that dark matter halos have been growing since equality. If t m j n > 2 t rec , 
then the computation is valid and we keep the result. 

For pa < there is an additional cutoff we must impose. One can see from Fig. |l| 
that halo virial density in such a universe continually decreases over time even as the 
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background matter density increases (following turnaround). Eventually we reach a 
point where halos virialize with densities less than that of the background, which is a 
signal of the breakdown of the spherical collapse model. 

The problem is that formally, the virial radius of the halo is larger than the radius 
of the "hole" that was cut out of the background and which is supposed to contain the 
"overdensity" . Thus, virialization really leads to a collision with the background, and 
not to a star forming isolated halo. Physically, this can be interpreted as a breakdown 
of the FRW approximation: sufficiently close to the crunch, the collapsing universe can 
no longer be sensibly described as homogeneous on average. 

Therefore, it makes sense to demand that the halo density be somewhat larger than 
the background density. To pick a number, we assume that under normal circumstances, 
a halo would settle into an isothermal configuration with a density proportional to 1/r 2 . 
Therefore the density on the edge of the halo will only be one third of the average 
density. If we demand a reasonable separation of scales between the density of the halo 
edge and that of the background matter, an additional factor of about 3 is not uncalled 
for. So we will impose the requirement that the average density of the halo be a factor 
of ten higher than the background matter density. 

We set the SFR to zero in any halos which do not satisfy this requirement. We 
implement this through a modification of the computation of t max - The true t max , 
computed at time t, is actually the minimum of the result above (Eqs. |2.37H2.38|) and 
the solution to 

Pvir(Uax) = 10p m (t) . (2.42) 



3. Results 

In this section, we apply our model to estimate the star formation rate, and the in- 
tegrated amount of star formation, in universes with various amounts of curvature, 
vacuum energy, and density contrast. 

3.1 Our universe 

We begin by considering our own universe. Figure |2| shows the SFR computed by our 
model as a function of time (left) and of redshift (right). At intermediate redshifts, 
our model is in good agreement with the data. At other redshifts, our model appears 
to be slightly off, but not by a large factor, considering its crudeness. For comparison, 
we show the SFR predicted by the models of HS [4], Nagamine et al. [15] and Hopkins 
and Beacom [16], which contain a larger number of fitting parameters. 

The shape of the SFR can be understood as a result of the interplay between 
competing forces. First, there will be no star formation until structure forms with 
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Figure 2: Comparison of our calculated SFR (solid) with the fossil model of Nagamine et 
al. [15] (dashed) and the HS model [4] (dotted). The right plot contains, in addition, the 
SFR of Hopkins and Beacom [16] (dash-dotted). The data points in the right plot are from 
astronomical observations and were compiled, along with references to their original sources, 



T vir > 10 4 K. Once a significant amount of structure begins to form at this critical 
temperature, however, the SFR rises rapidly and soon hits its peak. 

We can estimate the peak time, t pea k, of the SFR by asking when the typical mass 
scale that virializes is equal to M min : 



(There are ambiguities at the level of factors of order one; the factor a/2 was chosen 
to improve agreement with the true peak time of the SFR.) For the parameters of our 
universe, this calculation yields t pes ^ = 1.7 Gyr, close to the peak predicted by our 
model at about 1.6 Gyr. 

After t « tpeak, the SFR falls at least as fast as 1/t because t grav oc t. Once the 
typical virialized mass exceeds the cooling limit, M max , the SFR falls faster still. These 
considerations generalize to universes besides our own after accounting for the peculiar 
effects of each of the cosmological parameters. In particular, as explained below, we 
should be able to accurately estimate the peak time by this method for any universe 
where the effects of curvature or the cosmological constant are relatively weak. 



in [15] 




1 . 



(3.1) 
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Figure 3: The SFR for universes with greater cosmological constant (top) and more curvature 
(bottom) than ours. All other parameters remain fixed to the observed values, with curvature 
assumed to be open and just below the observational bound. 

3.2 Varying single parameters 

Next, we present several examples of star formation rates computed from our model 
for universes that differ from ours through a single parameter. 

Curvature and positive cosmological constant In Fig. |3], we show the SFR for 
our universe together with the SFRs computed for larger values of the cosmological 
constant (top) and larger amounts of curvature (bottom). One can see that increasing 
either of these two quantities results in a lowered star formation rate. 3 Essentially 

3 This differs from the conclusions obtained by Cline et al. [5] from an extrapolation of the Hernquist- 
Springel analytical model [4]. This model approximates all haloes in existence at time t as freshly 
formed. Moreover, it does not take into account the decreased relative overdensity of virialized halos in 
a vacuum-dominated cosmology. These approximations are good at positive redshift in our universe, 
but they do not remain valid after structure formation ceases. 
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Figure 4: The total stellar mass produced per comoving volume, as a function of open spatial 
curvature and the cosmological constant, (a) Curvature: Integrated SFR as a function of the 
number of efolds of inflation, relative to the observational bound on open spatial curvature. 
For AN > —2, star formation is unaffected. More curvature (less inflation) suppresses star 
formation, as expected. (The small irregularities in the asymptotic regime are numerical 
artifacts.) (b) Positive cosmological constant: p\ is given in Planck units; the observed value, 
PA ~ 1-25 x 10 -123 , is indicated by the position of the vertical axis. Decreasing p\ has no 
effect on star formation. Larger values interfere with structure formation, and star formation 
becomes suppressed, (c) Negative cosmological constant: Enhanced structure formation leads 
to more star formation for negative vacuum energy compared to positive vacuum energy of 
equal magnitude, (d) Positive and negative values of the cosmological constant shown on a 
linear plot. 

this is because structure formation as a whole is inhibited in these universes. Star 
formation is obstructed by the cosmological constant starting at the time t ~ = 
(3/A) 1 / 2 . Open spatial curvature suppresses star formation after the time of order 
t c = 0.533^ PeqG^eq) 3 , when fl = 0.5 in an open universe with p A = 0. Not shown 
are SFRs from universes which have a smaller cosmological constant than our universe, 
nor universes that are more spatially flat than the observational bound (A A" > 0). 
The SFRs for those choices of parameters are indistinguishable from the SFR for our 
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universe. 



In Fig. 4(a) , the integrated star formation rate is shown as a function of AiV. It is 
apparent that extra flatness beyond the current experimental bound does not change 
star formation. Indeed, the curve remains flat down to about AiV pa —2, showing that 
the universe could have been quite a bit more curved without affecting star formation. 
Integrated star formation is suppressed by a factor 2 for AiV = —3.5 and by a factor 
10 for AiV = -4. 

These numbers differ somewhat from the catastrophic boundary used in Ref. [17], 
AiV = —2.5. This is because the observational upper bound on flk has tightened from 
the value 20 x 10~ 3 used in Ref. [17] to 1.1 x 10~ 3 . The observationally allowed region 
corresponds to AiV > 0. Future probes of curvature are limited by cosmic variance to 
a sensitivity of $7& > 10~ 4 or worse, corresponding to AiV < 1 at most. The window 
for discovery of open curvature has narrowed, but it has not closed. 

The string landscape predicts a non-negligible probability for a future detection 
of curvature in this window. (This prediction is contingent on cosmological measures 
which do not reward volume expansion factors during slow-roll inflation, such as the 
causal diamond measure [2] or the scale factor measure [8]; such measures are favored 
for a number of reasons.) Depending on assumptions about the underlying probability 
distribution for AiV in the landscape (and on the details of the measure), this proba- 
bility may be less than 1% or as much as 10% [18]. However, the fact that curvature is 
already known to be much less than the anthropic bound makes it implausible to posit 
a strong landscape preference for small AiV, making a future detection rather unlikely. 

In Fig. [4(b)| , the integrated star formation rate is shown as a function of (positive) 
Pa- The observed value of p\ is right on the edge of the downward slope, where vacuum 
energy is beginning to affect star formation. Integrated star formation is suppressed 
by a factor 2 for p\ = 36pA and by a factor 10 for p\ = 225pA . To obtain a proba- 
bility distribution for p^ from this result, one needs to combine it with a cosmological 
measure [18]. 

Negative cosmological constant In Fig. [5] we see the SFR in some universes with 
negative cosmological constant. It is instructive to compare universes that differ only 
through the sign of the cosmological constant. The universe with negative vacuum 
energy will eventually reach a Big Crunch, where star formation ends along with ev- 
erything else. Yet, it will generally have a greater star formation rate than its partner 
with positive vacuum energy. This is because in the positive case, structure formation 
is severely hindered as soon as vacuum energy comes to dominate, at t pa t\, due to 
the exponential growth of the scale factor. In the negative case, this timescale roughly 
corresponds to the turnaround. Structure formation not only continues during this 
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(a) | p A | = pA 



(b) |p A | = 100p Ao 




(c) A7V = -4, |p A | =pa 

Figure 5: The SFR with negative cosmological constant. In each figure we consider two 
different universes whose parameters are identical except for the sign of p\. (a) Changing 
the sign of the observed cosmological constant would have made little difference to the SFR. 
However, in universes with a higher magnitude of cosmological constant (b), or more highly 
curved universes (c), a negative cosmological constant allows more star formation than a 
positive one of equal magnitude. The apparent symmetry in (c) is accidental. 



phase and during the early collapse phase, but is enhanced by the condensation of the 
background. Star formation is eventually cut off by our requirement that the virial 
density exceed the background density by a sufficient factor, Eq. Q2.42Q . This is the 
origin of the precipitous decline in Fig. [|(b) and (c). 

The integrated star formation rate as a function of negative values of the cosmo- 
logical constant is shown in Fig. [|. As in the positive case, sufficiently small values of 
| Pa | do not affect star formation. A universe with large negative pa will crunch before 
stars can form. For intermediate magnitudes of pa, there is more star formation in 
negative than in positive cosmological constant universes. The amount of this excess 
is quite sensitive to our cutoff prescription near the crunch, Eq. ( ggp . With a more 
lenient cutoff, there would be even more star formation in the final collapse phase. A 
stronger cutoff, such as the condition that p vu > 100 p m , would eliminate the excess 
entirely. Clearly, a more detailed study of structure and star formation in collapsing 
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Figure 6: The SFR for different values of the primordial density contrast, Q. The upper 
plot shows the great sensitivity of the SFR to relatively small changes in Q. Compton cooling 
plays no role here. The lower plot shows a larger range of values. The universal slope of 
— 1 in the log- log plot corresponds to the 1/t drop-off in the Compton cooling regime. When 
Compton cooling becomes ineffective at around 1 Gyr, most of the SFRs experience a sudden 
drop. The spike at early times for large values of Q arises from the sudden influx of baryons 
after recombination. 



universes would be desirable. 

Density contrast In Fig. [] we see the effect of increasing or decreasing the amplitude 
of primordial density perturbations, Q. Note that even relatively small changes in Q 
have a drastic impact on the SFR. Increasing Q has the effect of accelerating structure 
formation, thus resulting in more halos forming earlier. Earlier formation times mean 
higher densities, which in turn means higher star formation rates. 

Of course, higher Q also leads to larger halo masses. One might expect that cooling 
becomes an issue in these high-mass halos, and indeed for high values of Q the SFR 
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Figure 7: In our model, total star formation increases with the density contrast. This 
increase continues until we can no longer trust our SFR due to the dominance of large-scale 
structure at recombination. As discussed in the text, a more refined model may find star 
formation inhibited in parts of the range shown here. 

drops to zero once Compton cooling fails, at t = r comp . However, for t < r comp , cooling 
does not limit star formation, and so in high Q universes there tends to be a very large 
amount of star formation at early times. For the highest values of Q that we consider, 
there is a spike in the star formation rate at t = 2t vec coming from the sudden infall 
of baryons to the dark matter halos. The spike is short-lived and disappears once that 
initial supply of gas is used up. 

In Fig. we show the integral of the SFR as a function of Q. The increase in 
total star formation is roughly logarithmic, and in our model the increase keeps going 
until most structure would form right at recombination, which signals a major regime 
change and a likely breakdown of our model. 

A universe with high Q is a very different place from the one we live in: the densities 
and masses of galaxies are both extremely high. It is possible that qualitatively new 
phenomena, which we are not modelling, suppress star formation in this regime. For 
example, no detailed study of cooling flows, fragmentation, and star formation in the 
Compton cooling regime has been performed. We have also neglected the effects of 
enhanced black hole formation. 

3.3 Varying multiple parameters 

We can also vary multiple parameters at once. For example, let us hold only Q fixed 
and vary both curvature and (positive) vacuum energy. The integrated star formation 
is is shown in Fig. ||. If either curvature or pa get large, then structure formation 
is suppressed and star formation goes down. When the universe has both small pa 
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Figure 8: The integrated SFR as a function of curvature (AN) and positive p\, shown in 
a 3D plot (a) and in a contour plot (b). (The strength of density perturbations, Q, is fixed 
to the observed value.) Analytically one expects star formation to be suppressed if t pea k > t c 
o r ^peak > where i pea k is the peak star formation time in the absence of vacuum energy 
and curvature. These catastrophic boundaries, log 10 pa = —121 and AiV = —3.5, are shown 
in (b). 




Figure 9: The integrated SFR as a function of positive p\ > and Q, at fixed spatial 
curvature, AiV = . The contour plot (b) also shows the catastrophic boundary i pea k(Q) > *A 
expected analytically. 



and small curvature, then structure formation is not suppressed and the SFR is nearly 
identical to that of our universe. Here "large" and "small" should be understood in 
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terms of the relation between t\, t c , and t pea k computed in a flat, pa = universe 
(an unsuppressed universe). In the case we are considering, namely Q = Q , the 
unsuppressed peak time is t pea k = 1-7 Gyr. The conditions t peak = t\ and t pea k = 
translate into pa ~ 100 pa and AN w —3.5 . These lines are marked in Fig. |8(b)| , and 
one sees that this is a good approximation to the boundary of the star-forming region. 
(We will see below that this approximation does not continue to hold if both curvature 
and Q are allowed to vary.) 





(a) 



(b) 



Figure 10: The integrated SFR as a function of AN and Q, at fixed p\ = p\ Q . The 
contour plot (b) also shows the boundary i pea k(Q) > t c expected analytically, but this analytic 
boundary is not a good approximation to the true catastrophic boundary. 



If we also start to vary Q, the story gets only slightly more complicated. Increasing 
Q causes structure formation to happen earlier, resulting in more star formation. In the 
multivariate picture, then, it is helpful to think of Q as changing the unsuppressed t pea k- 
For the unsuppressed case, the universe can be approximated as matter-dominated. 
Looking at Eq. |3.1| and recalling that G ~ a oc 

t 2/3 

in a matter-dominated universe, 
and neglecting for the sake of computational ease the dependence on M min , we see that 
£ pea k oc Q~ 3 ^ 2 - Now we can quantitatively assess the dependence of star formation on 
Q and pa, for instance. We know that £a = £ pea k when Q = Q and pa = 100 pa , 
and this condition should generally mark the boundary between star formation and no 
star formation. Therefore we can deduce that the boundary in the two-dimensional 
parameter space is the line log 10 pA/pA = 2 + 31og 10 Q/Q ■ This line is shown in 
Fig. P(b) , where it provides a good approximation to the failure of star formation. 

We can repeat the same analysis for the case of fixed cosmological constant with 
varying curvature and perturbation amplitude. In Fig. ITTJ we show the integrated star 
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Figure 11: The linear growth factor for three different universes: a flat p\ 7^ universe 
(solid), a curved p\ = universe (dot-dashed), and a flat, p\ = (unsuppressed) universe 
(dshed). The former two universes have been chosen so that t c = t\ = i cr i t , and the time 
axis is given in units of i C rit- In the p\ case, the growth factor can be well-approximated 
as being equal to the unsuppressed case up until t = t cr it, and then abruptly becoming flat. 
This approximation is manifestly worse for the curvature case. In the curvature case, the 
growth factor starts to deviate from the unsuppressed case well before t = t c , and continues 
to grow well after t = t c . This effect makes the case of strong curvature harder to approximate 
analytically. 



formation rate as a function of AiV and Q for fixed pa = Pa - The dashed line in 
Fig. |10(b)| marks the boundary t c = t pea k, where the time of curvature domination 
equals the time of the unsuppressed peak (again computed according to Eq. |3.1| ). In 
this case the line is given by (2 AiV + 7) log 10 e = — log 10 Q/Qo, owing to the fact that 
t c oc exp(3AiV). Unlike the case of pa, this line obviously does not mark the transition 
from structure formation to no structure formation. The reason has to do with the 
details of the nature of the suppression coming from curvature versus that coming from 
the cosmological constant. In the case of the cosmological constant, it is a very good 
approximation to say that structure formation proceeds as if there were no cosmological 
constant up until t = t\, whereupon it stops suddenly. Curvature is far more gradual: 
its effects begin well before t c and structure formation continues well after t c . Fig. |H] 
illustrates this point clearly. There we see that while perturbations are suppressed 
in both the high curvature and large pa cases, it is only in the pa case that the naive 
approximation of a "sudden end" to structure formation is valid. For the curvature case 
there is no simple approximation, and the full model is necessary to give an accurate 
assessment of the situation. 
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